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ABSTRACT 

We examine the recoil velocity induced by the superposition of the magnetic dipole 
and quadrupole radiation from a pulsar/magnetar born with rapid rotation. The resul- 
tant velocity depends on not the magnitude, but rather the ratio of the two moments 
and their geometrical configuration. The model does not necessarily lead to high spatial 
velocity for a magnetar with a strong magnetic field, which is consistent with the recent 
observational upper bound. The maximum velocity predicted with this model is slightly 
smaller than that of observed fast-moving pulsars. 

Subject headings: stars: magnetic field — stars: neutron — pulsars: general 
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1. INTRODUCTION 

The surface magnetic field strength B s of a pulsar is conventionally estimated by matching the 
rotational energy loss rate with the mag netic dipole radiation rate, that is, B s « (Zc'IPP) 1 / 2 / '(2 v / 2vr^) ! 
where I is the inertial moment, R s is the stellar radius, P is the spin period, and P is the time deriva- 
tive of the spin period. The precision of this approximation is only at the order of magnitude level 
because actual energy loss is not well described by magnetic dipole radiation in a vacuum. A more 
realistic model with current flows and radiation losses is required, but has not yet been established. 
A simple estimate provides B s ~ 10 12 G for typical radio and X-ray pulsars, and B s ~ 10 13 -10 15 G 
for magnetars, although the level of the approximation must be noted. Dynamo action in a rapidly 
rotating proto-neutron star with P 1 ms is propose d as a mechanism for this amplifica t ion b y 
2-3 orders of magnitude (see e . g.. iThompson fc Duncan! (|1993l ); lBonanno. Rezzolla k, Urpinl (|2003l ); 
Bonanno. Urpin &: Belvedere! (|2006l )). Actual upper limit of B s generated in the convective proto- 
neutron star is estimated as 10 15 — 10 16 G, beyond w hich all sorts of instabilities are suppressed by 



strong magnetic fields (jMirailes. Pons &: Urpinl I2002I ). 



Recent numerical simulations of dynamo action can be used to study the large-scale fields in 
fully convective r otating stars. For examp le, non-axisymmetric fields are generated in the case of 
uniform rotation (jChabrier &: Kiikerll2006l ). while mostly axisymmetric fields with a mixture of the 
first f ew multipoles are formed in the case of a differentially rotating star ljDobler. Stix Brandenburg 
2006). The results may not directly apply to pulsars or magnetars, but suggest that the magnetic 
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field configuration of neutron stars may not be an ordered dipole. If there are higher-order multi- 
poles, these will also contr ibute to the ra diation loss. The upper bounds on their surface magnetic 
fields are rather loose. See iKrolik I (|199ll ) for a discussion of the magnetic fields of millisecond pul- 
sars. The magnetic field strength B[ m relevant to the multipole moment of order (I, m) is limited to 
Bim < B s /(mR s tt/c) 1 - 1 , where O = 2vr/P is angular velocity, and the radiation of each multipole 
Lim ~ c(mR s Q I c) 2l+2 (Bi m R s ) 2 is assumed to be smaller than that of a dipole. Thus, a model 
with complex magnetic configuration at surface Bi m > B s (I > 1) is allowed because of the small 
factor R s £l/c <C 1 for observed stars. 

Some proto-neutron stars are conjectured to be born in hypothetical extreme state of rapid 
rotation P ~ 1 ms with an ultra strong magnetic field B s ~ 10 15 G. Is there any remaining evi- 
dence of this stage? The proper motion can possibly be used as a probe. Several kick mechanisms 
operative at the core b ounce of a supernova explosion have been prop osed to date: anisotropic emis- 



sions of neutrinos ( e.g. jArras Lail (jl999l ); lFrver fe K uscnko (2006)), hydrodynamical waves (e.g., 



Scheck et al.l (|2006l )). and MHD effects (e.g.. ISawai. Kotake fc Yamadal (120081 )). These mechanisms 
operate on a dynamical timescale of the order of milliseconds or the cooling timescale of ~ 10 s. If 
the strong magnetic fields are generated on a longer timescale, some natal kick mechanisms involved 
the magnetic-field-driven anisotropy do not work effectively. Recoil driven by electromagnetic radi- 
ation, which is operative on a longer spin down timescale of ^ 10 3 (P/1 15 G)~ 2 (P ? ;/lms) 2 s, has been 
propo sed as a post-natal kick mechanism ( Harrison fe Tademarulll975 ) (see also Lai. Chernoff h Cordes 
(|200ll ) for the corrected expression). In their model, an oblique dipole moment displaced by a dis- 
tance s from the stellar center rotates. This causes the radiation of higher order multipoles, whose 
superposition is generally asymmetric in the spin direction, leading to the kick velocity. In the off- 
center model, the quadrupole field B2 of order B2 ~ (s/R s ) x B\ ~ B\ is involved. It is interesting 
to study the case where B2 S> B\, because the constraint of the higher order component by the 
radiation is very weak, for example, B2 ~ (c/(R s £l)) x B\ 3> B\. In this paper, we revisit the kick 
velocity induced by electromagnetic radiation from a magnetized rotating star with both dipole 
and quadrupole fields, in which a larger quadrupole field B2 3> B\ at the surface is allowed. 

This paper is organized as follows. In Section 2, we present the radiation from rotating dipole 
and quadrupole moments in vacuum. Their field strength and inclination angle with respect to the 
spin axis are arbitrary. We also compare our model with the off-centered dipole model. In Section 
3, we evaluate the maximum kick velocity as a recoil of momentum radiation. Section 4 presents 
our conclusions. 



2. MODEL 



2.1. Electromagnetic Fields 



We consider electromagnetic fields outside a rotating object with angular frequency f2; the 
object has a magnetic dipole and quadrupole moments. The dipole moment is denoted by fi, and 
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the direction is inclined from the spin axis by xi- Quadrupole moment is denoted by Q and the 
inclination angle of the symmetric axis is X2 from the spin axis. The electromagnetic fields outside 
the rotating magnetized object are desc ribed by the magnetic mupltipoles of order I = 1, 2, \m\ < I, 
for which E r = E • r = (j Jackson Ill975l ). The explicit components produced by the rotating dipole 
moment are given by 
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where \ = <p — VL(t — r / c). The function P r f{x) is the associated Legendre function and the prime 
denotes the derivative with respect to x. The function S n , which is derived from the spherical 
Hankel function, is a polynomial of £ = £lr/c, and is explicitly written as 



Si (0 = i-*£, 



(6) 
(7) 



We here use convenient complex expressions in eqs. (PQ)— dSD an d the actual fields are the real part. 
Near the stellar surface R s < r <C c/f2, magnetic field for eqs.(HD-([3]) at the phase e lX = 1 reduces 
to 
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It is clear that the field near the origin represents a magnetic dipole inclined by the angle xi> which 
rotates in the azimuthal direction with (f> = £lt. 

The electromagnetic fields for a rotating magnetic quadrupole are similarly described by 
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where 

s 3 (o = i-»e- g£ 2 , ( 14 ) 

5 4 (0 = l-^-k 2 + ^ 3 . (15) 

The phase in eqs. (|9j)- (ll3p is shifted by A2 = A + 5, because the meridian plane in which the 
symmetric axis of the quadrupole is located may differ by the azimuthal angle 5 from that of the 
dipole. The near-field of eqs.([9l)- (fTTj) for R s < r -C c/fl at the phase e* A2 = 1 is 

B = ~ X2)e r - ^P%{9 - X2)e e = -V (^P 2 > - »)) • (16) 

Thus, the magnetic fields given by eqs. (|9|)- (|lip are those of a rotating quadrupole, whose inclination 
angle is X2- 

We compare the combination of dipole and quadrupole fields with the case of a pure dipole. 
A snapshot of almost-closed magnetic field lines near the light cylinder is shown in Fig. [U Both 
inclination angles are the same xi = X2 = t/4, but the meridian planes are perpendicular, that is, 
5 = 7r/2. Field strength is set as Q = 0.2/ic/O. It is clear that the quadrupole field is added to 
the dipole one. The quadrupole field increases more rapidly with the decrease of the radius r, and 
dominates for r < r q ~ 0.2c/fi for the model parameter, since B\ ~ fi/r s and B2 ~ Q/r 4 . 




Fig. 1. — Closed magnetic field lines in the cases of pure dipole (left) and dipole plus quadrupole 
(right). The magnetic axis is inclined by angle xi = X2 = from the spin axis z, and the 
azimuthal angle between moments is 5 = n/2. A sphere of radius 0.2c/O at the origin is also 
shown. Distance is scaled by c/Q. 
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2.2. Radiation 

The radiation energy per unit time is obtained by integrating the time-averaged Poynting flux 
over the solid angle at the wave zone r ^> c/Cl. The luminosity for a combination of electromagnetic 
fields described by eqs. dH)-© and eqs. ([9j)- (fT3|) is given by 

L = J — (E x B) • e r r 2 sin0d0# = sin 2 X i + sin 2 2 X2 + sin 4 X2 . (17) 

The luminosity is the sum of the contributions from multipole radiation. Our model consists 
of three components, the magnetic dipole radiation Mi 1 specified by spherical harmonics index 
(l,m) = (1, 1), and the quadrupole radiation il^i and M^i- They correspond to the first, second 
and third terms in eq. (1171) , The third term is larger than the second term roughly by a factor 
m 6 = 2 6 , which comes from the frequency of time variation. 

The linear momentum radiated per unit time in the direction z is similarly calculated as 

/I u0f2 5 
— (E x B) • e z r 2 sin 9d6d<f> = ^ Q 5 sin xi sin 2x2 sin 5. (18) 

The net flux arises from the interference of two multipoles, namely, the magnetic dipole M\ t \ 
and the quadrupole M21. The angle xi governs the radiation strength of each multiple /, while 
the angle S governs the interference. The most efficient configuration is realized when the two 
magnetic multipole moments are orthogonal, 5 = ir/2. On the other hand, when both of the 
multipole moments lie in the same meridian plane (i.e., 5 = 0), the net linear momentum vanishes. 
This property can be understood from the fact that radiative electromagnetic fields in vacuum 
are expressed by the spherical Hankel function hi and the asymptotic form for £ = Or/c 3> 1 is 
hi ~ exp[i(£— lir/2)]/r for the multipole I. There is a phase shift tt/2 between dipole and quadrupole 
fields, and this shift is important in the wave interference. 



2.3. Comparison 



We compare our result with the off-center dipole model ( Harrison &: Tademarulll975 : Lai. Chernoff Cordes 



200ll ) . The rates of energy and linear momentum are written in term of the magnetic dipole moment 



^4>i Mz) m cylindrical coordinate and distance s from the spin axis as follows: 

L = l? (/4+/4) + ^ 2 /^- (19) 

The first term is the magnetic dipole radiation M\ t \. Correspondence to our expression is clear by 
replacing /i^ + jU^ = /U 2 sin 2 xi- The second term is derived from the sum of electric dipole radiation 
Ei 1 and magnetic quadrupole radiation M21. Their contributions are 6 s 2 /i 2 /(6c 5 ) by Ei \ and 
s 2 fi 2 1 (Wc 5 ) by Mjj,!, respectively. The parameter in the off-center dipole model corresponds 
to Q sin 2x2 = 4s/^ except for a complex phase factor. There is a constraint on the quadrupole 
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moment Q as Qsin2x2 < 4//-R s cosxi> since s < R s . In our model, it is possible to consider the 
case of Q S> ^R s in magnitude. 



The linear momentum in the off-center dipole model is evaluated as (|Lai. Chernoff Cordes 

200ll ) 

F = 8 "W. (20) 

loc 

Net linear momentum flux arises from two types of interference. One is between magnetic dipole 
radiation and electric dipole radiation E\\. The other is between magnetic dipole radiation 
M\ t i and magnetic quadrupole radiation M 2t i- These contributions are expressed by il°sfj,^fj, z / (3c 5 ) 
and f2 5 s/x<^ 2 /(5c 5 ), respectively. The latter reduces to eq. ([T5]) if sfi z = Q sin 2%2/4 and /i^ = 
/isinxi sin S. 

Although there is a slight difference in the radiative components between the off-center dipole 
and dipole-quadrupole models, both formulae for eqs. (fT7"|) . (fl~8)) and eqs. (fT9j) . ([20]) are parameter- 
ized as 

L = a^ + /3^-, (21) 

F "10^5-' (22) 

where a, (3 and 7 are dimensionless numbers that depend on only the geometrical configuration. 
The typical values are listed in Table [1] for the simple assumption that sinxz,sin<5 — > l/\/2, that is, 
the directional average of (sin 2 xi) = (sin 2 S) = 1/2. It is clear that the coefficient /3 in our model 
is considerably larger than that in the off-center model. This comes from the radiation of m = 2. 



3. EVOLUTION 

We next consider the evolution of spin and kinetic velocity. The angular velocity Q(t) is 
determined by equating the loss rate of rotational energy with the luminosity L in eq. (f2Tj) . and 
the velocity V(t) is determined from the momentum emission F in eq. (|22p . In terms of the mass 
M and inertial moment /, we have 

m O = -a^--/3^-, (23) 

cr c° 



Table 1: Comparison of models. 



Model 


Multipole 


a 


P 


7 


7 /(a/3) 1 /2 


Off-center dipole 


M 1;1 ,M 2:1 ,E 1:1 


0.33 


0.83 xl0~ 2 


0.47 


9.0 


Dipole-quadrupole 


M lil ,M 2 ,i,M 2 ,2 


0.33 


0.10 


0.18 


0.97 
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MV = - 1 L^^. (24) 



7 mQ^ 5 
10 c 5 

By using the approximation / = 2Mi? 2 /5, where R s is the stellar radius, the magnitude of the 
velocity gained from the initial angular velocity Oj is given by 



AV = £ at+KQiriw * 1 - AV - s (n,R - f x ' 2 [x - ^ x] ' <25> 

where X = (J3/a)^' 2 QSli/ \i and the present angular velocity f^o = is used in the last inequality. 
The function AV* is determined by the ratio Q/(fiR s ) of the two multipole moments for the 
fixed geometrical configuration and the initial angular velocity f2j. Two limiting cases of AV* are 
approximated as 



AV* I rk{A)(^) 3 forO<X«l 



7 _Q_\ l ^iRs 

25/3 [fiR, 



(26) 
for X > 1. 



The value AV* increases as the ratio Q/(fiR s ) increases, while Q/(fiR s ) <C 1, but begins to 
decrease for Q/(fiR s ) — > oo. Thus, it has a maximum with respect to the magnetic moment ratio: 

^ K9 . 2xl „-3 7 (Ml)' at « M L5 ( "V'VfMiV 1 . (27) 



The magnetic moment ratio at the maximum means that the quadrupole field B2 ~ Q/Rs is 
stronger than the dipole field B\ ~ n/Rg at the surface. The energy loss rate L; of each multipole 
is approximately the same at the beginning, L 2 = (3Q 2 nf/c 5 « 2.3 x a^ 2 ^/ c 3 = 2 .3ii, but the 
contribution of Li and becomes smaller as O is decreased. The velocity using the canonical values 
is evaluated as 2 

AV* « 120 f/U x —J— km s- 1 (28) 

For the off-center dipole model, V* ~ 10 3 (Pj/lms)~ 2 km s -1 is allowed for an initially rapid 
rotator, the initial period Pj = 1 ms, using typical values given in Table [TJ On the other hand, 
the typical value is small, V* ~ 10 2 (p/lms) -2 km s _1 , for the dipole-quadrupole model. The 
difference comes from the presence of radiation of m = 2, which causes efficient energy loss, as 
discussed in the previous section. Nevertheless, extremely high velocity is possible for a specific 
configuration even in the present model. Small (3 corresponds to high velocity. For small X2 in eq. 
(fl7l) . we have /3 = sin 2 2X2/I6O. Because a = 2sin 2 xi/3, 7 = sinxi sin 2x2/2 for sin 5 = 1; the 
combination of parameters reduces to 7/(a/3) 1 / 2 = 7.7. The resultant kick velocity increases up to 
~ 930(Pj/lms) -2 km s _1 . This optimal case corresponds to the magnetic configuration with an 
inclined dipole and a nearly axially symmetric quadrupole. The ratio of the moments is Q/{jj,R s ) ~ 
74(sin xi/ sin X2)(P/lms). 

Time evolution of spin and velocity is calculated for the optimized relation ()27|) . Once the 
quadrupole field strength is fixed, the evolution of Q(t) in eq. (|23l) is scaled by characteristic time 
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i* of the dipole radiation loss for the initial angular velocity f2 

Ic 3 f a \ -l / Pi 



2*/ Pc 



0.8 



a 
033 



(. 



lms J \10 31 Gcm 2 



-2 



yr , 



(29) 



where magnetic dipole field at the surface is chosen as B\ ~ 10 G. Figure [2] shows the evolution of 
as a function of r = i/i*. The ratio of quadrupole to the total energy loss rate, L2/(L\ +L 2 ) 
is also plotted. The ratio at t = is approximately 0.7 because of L2 ~ 2.3Li, but monotonically 
decreases. At t = i*, the angular velocity becomes Q ~ 0.5ilj and the contribution of quadrupole 
radiation also decreases as L 2 ~ 0.5Li. The velocity V(t) normalized by the terminal one (eq.(|28|)) 
is also shown in Figj2j The magnitude attains to almost terminal value, V ~ 0.8AI4 before t = t*. 
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Fig. 2. — Time evolution of velocity, spin and energy loss rates as a function of dimensionless 
time t = t/t*. Solid line represents the linear velocity normalized by the terminal one A V*, and 
dashed line the angular velocity normalized by initial one Oj. Dotted line is the radiation loss ratio 
L 2 /(L X + L 2 ). 



4. CONCLUSIONS 



Magnetic field strength itself is critical in most kick mechanis ms. For example, B s > 10 15 
G at the surface is required in asymmetric neu trino emission (e.g., Arras &: Lail (|1999l )). as well 



as in asymmetric magnetized core collapse (e.g., ISawai. Kotake &; Yamadal (|2008l )). The resultant 
velocity increases with the field strength because the asymmetry arises from the magnetic field. 
Magnetars are therefore expected to have high velocity if one of these mechanisms is operative. 
Recent observations do not support the high velocity. Rather, the upper limit of the transverse 
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velocity v± has been reported, altho ugh there is uncert ainty in the value. For example, v± ~ 
210 km s" 1 for AXP XTEJ1810-197(lHelfand et al.l 120071^ . v ± < 1300 km s" 1 for SGR 1900+14 



(Ka plan et al.ll2009l : be Luca et al.ll2009h and v± < 930 km s" 1 for A XP lE2259+586(IKaplan et al 



2009). For fast moving pulsars, PSR2224 +45 (v \ > 800 km s 1 ICordes. Romani &: Lundgren 



(jl993j)) and B1508+55 (v ± ~ 1000 km s" 1 IChatteriee et all (|2005j)) have been reported. These 
magnetic fields are quite ordinary, B s = 2.7 x 10 12 G and 2.0 x 10 12 G, respectively. Thus, there is 
no clear correlation between the field strength and the velocity in the present sample. 



The electromagnetic rocket mechanism considered in lHarrison &: Tademarul (j 19751 ) and in this 
paper does not depend on field strength if the spin evolution is determined from the radiation loss. 
In our model, the ratio of dipole and quadrupole moments is important. The condition for high 
velocity is that the quadrupole field is large enough in magnitude for the radiation loss to be of the 
same order as the dipole field. The velocity also depends on the geometrical configuration of the 
multipole moments, that is, each inclination angle from the spin axis and the angle between the axes 
of symmetry of the moment. Assuming that the directions of moments are random, and that they 
are equally likely to be oriented in any direction, it is found that the mean velocity with respect 
to the configuration is not so large, ~ 120(Pj/lms) -2 km s , for the optimized dipole-quadrupole 
ratio. The maximum velocity is realized for a specific configuration in which the inclination angle 
of the quadrupole moment is small, and the meridian plane in which the quadrupole moment lies 
is perpendicular to the plane of the dipole. The velocity increases up to ~ 930(Pj/lms) -2 km s . 
This value is slightly smaller than the maximum observed velocity of a pulsar. 

The configuration is unknown, and is closely rela t ed to the origin of the magnetic field, dynamo 
or fossil. Nevertheless, iBonanno. Urpin Belvedere! (|2006l ) reported interesting results within the 
mean-field dynamo theory. They argued that strong large-scale and weak small-scale fields are 
generated only in a star with a very short initial period, that is, the Rossy number is small, and 
that the maximum strength decreases and small-scale fields become dominant with the decrease 
of the initial period. Thus, magnetars may have an ordered dipole with a strong field, while some 
pulsars may have rather irregular fields with higher multipoles. Through the superposition of higher 
multipoles, pulsars in general come to have a larger radiation recoil velocity than magnetars. 

Finally, if the kick velocity of pulsars and magnetars is governed by the same mechanism, it 
either should not simply depend on magnetic field, or should depend on only the configuration. 
The latter possibility was explored here. Present argument is recognized as the order of magnitude 
level due to the rotating model in vacuum. Further improvement of the magnetosphere will be of 
importance to explore the idea. 
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